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ABSTRACT 

Previous analyses of magnetospheric accretion and outflow in classical T Tauri 
stars (CTTSs), within the context of both the X-wind model and other theoretical 
scenarios, have assumed a dipolar geometry for the stellar magnetic field if it 
were not perturbed by the presence of an accreting, electrically conducting disk. 
However, CTTS surveys reveal that accretion hot spots cover a small fraction 
of the stellar surface, and that the net field polarization on the stellar surface 
is small. Both facts imply that the magnetic field generated by the star has 
a complex non-dipolar structure. To address this discrepancy between theory 
and observations, we re-examine X-wind theory without the dipole constraint. 
Using simple physical arguments based on the concept of trapped flux, we show 
that a dipole configuration is in fact not essential. Independent of the precise 
geometry of the stellar magnetosphere, the requirement for a certain level of 
trapped flux predicts a definite relationship among various CTTS observables. 
Moreover, superposition of multipole stellar fields naturally yield small observed 
hot-spot covering fractions and small net surface polarizations. The generalized 
X-wind picture remains viable under these conditions, with the outflow from a 
small annulus near the inner disk edge little affected by the modified geometry, 
but with inflow highly dependent on the details of how the emergent stellar flux 
is linked and trapped by the inner disk regions. Our model is consistent with 
data, including recent spectropolarimetric measurements of the hot spot sizes 
and field strengths in V2129 Oph and BP Tau. 

Subject headings: stars: low mass - stars: pre-main-sequence -circumstellar mat- 
ter - stars: formation 
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1. Introduction 

Various theoretical models have been proposed for the physical mechanisms driving the 
accretion and outflow processes in classical T Tauri stars (CTTSs), with the most popular 
being perhaps Blandford and Payne's (1982) pioneering study of magnetocentrifugally driven 
winds from Keplerian disks. However, the X-wind model has gained credence in recent years 
for a variety of theoretical and observational reasons (see also the discussion in §5). It is 
instructive at the outset to summarize and compare the different CTTS models, to review 
the rationale for X-wind theory and the evidence in its favor, and to motivate our subsequent 
generalization of this picture (see also Shu et al. 2000). 

The first accretion models of CTTSs proposed that they were surrounded by Keplerian 
disks, which extend all the way to the stellar surface, with accretion occuring through an 
equatorial boundary layer (e.g., Bertout 1987; Kenyon & Hartmann 1987; Bertout et al. 
1988). Motivated by the observational finding that the driver behind some well-known 
bipolar outflow sources are neutral winds containing H I and CO (Lizano et al. 1988), a 
combination not seen in the interstellar medium but present in the photospheres of cool 
stars, Shu et al. (1988) suggested that the driving outflows are caused by boundary-layer 
disk-accretion onto a strongly magnetized young stellar object (YSO). The accretion provides 
a rationale for why a protostar might rotate near break-up, and Hartmann & MacGregor 
(1982) had already shown that magnetized stars rotating near break-up could shed matter 
and angular momentum extremely efficiently through their equatorial zones. However, Shu et 
al. (1988) obtained preliminary indications that streamline collimation in this kind of model 
was absent, or extremely slow, and therefore they speculated stellar jets were associated 
with "ordinary" stellar winds confined to flow along the rotation axis by the more powerful 
"extraordinary" centrifugally-driven outflow that they later called an X-wind. 

When a disk abuts a fully convective star, Ekman pumping (internal circulation caused 
by slight pressure differences in a field of differential rotation) to the equatorial boundary 
layer adjoining the star and the disk eventually causes the entire star to spin near break-up 
if there are no countervailing spin-down torques (Galli 1990; see Fig. 10 of Shu et al. 1993). 
CTTSs are fully convective, yet they usually rotate at velocities one order of magnitude 
below break-up (e.g., Vogel & Kuhi 1981; Hartmann & Stauffer 1989; Bouvier et al. 1993, 
2007). Shu et al. (1988) therefore speculated that CTTSs are young stars in which the spin- 
down by ordinary stellar winds overcame the spin-up toques of the viscous boundary-layer, 
and that CTTSs would not have strongly collimated outflows. 

The latter expectation turned out to be false (see, e.g., Edwards et al. 1993). Never- 
theless, Matt & Pudritz (2005, 2008a, b) have recently resurrected the idea that magnetized 
stellar winds might provide sufficient torque to explain the slow rotation of T Tauri stars 
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and, perhaps, be responsible for a part of the observed optical jets from YSOs. The lat- 
ter proposal using coronal gas to launch the stellar winds has, however, foundered on the 
constraints provided by X-ray observations of T Tauri stars (Bisnovatyi-Kogan & Lamzin 
1977, DeCampli 1981, Matt & Pudritz 2008c), and the finding that CTTS stellar winds 
appear to be an order of magnitude cooler (Johns-Krull & Herczeg 2007) than suggested 
recently (Dupree et al. 2005). The issues, therefore, became (1) how to accrete high angular 
momentum disk material onto a CTTS while keeping the latter rotating slowly, (2) how si- 
multaneously to generate winds efficiently from the star+disk system, and (3) how to make 
the outflowing wind appear as a jet despite the slow streamline collimation? 

Konigl (1991) suggested a solution for the accretion part of this problem by adopting 
the theory developed by Ghosh & Lamb (1978; 1979a, b). A strong stellar magnetosphere is 
assumed to truncate the disk some distance from the stellar surface, with accreting material 
flowing onto the star not through a boundary layer but via magnetic field lines threading 
the disk. The angular momentum of the star is then regulated by the interaction between 
the stellar magnetosphere and the disk. Field lines originating in the star and with disk 
footpoints within the Keplerian corotation radius (Rx) are dragged forward by the disk gas, 
relative to the star, and thus tend to spin the star up (at the expense of the disk gas). Stellar 
field lines threading the slower rotating disk beyond Rx, on the other hand, tend to spin the 
star down. The slow rotation of CTTS then arise if the spin-down caused by the outer disk 
outweighs the spin-up by the inner parts. Large turbulent diffusivities must prevail to allow 
magnetic fields to slip through the gas in such a picture. 

In a series of papers with a precedent in the suggestions of Arons (1986) and Camenzind 
(1990), Shu and collaborators (Shu et al. 1994a, b; Najita & Shu 1994; Ostriker & Shu 1995, 
hereafter OS95; Shu et al. 1995) embraced the idea that strong YSO magnetospheres might 
truncate the disk before it abuts the stellar surface, but they pointed out a number of 
difficulties with the specific proposal of Konigl (1991). The resolution of these difficulties 
turned out to provide solutions for each of the issues (1) to (3) listed above. 

First, the picture painted by Ghosh & Lamb and Konigl changes considerably if turbu- 
lent resistivities are not large, but small, i.e., if field diffusion is competitive with advective 
flow not on dynamical time scales but secular ones. In such cases, the long-term processes 
of angular-monetum transport outwards and mass-transport inwards dominate over notions 
of "ram-pressure balance". Second, stellar fields strong enough to truncate an electrically 
conducting disk are automatically also strong enough to drive a magnetocentrifugal outflow 
along the outermost flux tubes of the stellar field, which are opened into an X-wind. The 
magnetic torques in the wind cause the outflowing matter to gain angular momentum at the 
expense of the material still connected to it by field lines threading through the disk. The 
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back reactions to the X-wind and funnel-flow give a pinch of the exterior field lines inward 
and the interior field lines outward towards a common mid-point Rx, with a net trapped flux, 
that gives the X-wind model its name. This mid-point is both in Keplerian rotation and in 
corotation with the star, i.e., Q(Rx) = (GM*/ R x ) 1 / 2 = fi*, a condition that came to be 
called disk locking. Third, isodensity contours become cylindrically stratified very quickly 
after the gas accelerates from the X-region (Najita & Shu 1994; Shu et al. 1995), yielding 
the optical illusion in the emission of forbidden-line and radio emission that X- winds achieve 
jet-like collimation close to the base of the flow (Shang et al. 2002, 2004). In actual prac- 
tice, streamline collimation is logarithmically slow, which has observable consequences for 
position- velocity diagrams obtained from long-slit spectrograms (e.g., Pyo et al. 2006). 

Many numerical simulations have also been devoted to the X-wind/funnel flow problem, 
the most successful being that of Romanova et al. (2008). Indeed, the progress made by 
the simulations is most concisely revealed by examining why they succeeded in obtaining 
the simultaneous existence of X- winds and funnel flows over extended durations when others 
failed. As stated in the previous paragraph and detailed in §2, the crucial concept in X- 
wind/funnel flow theory is that of the trapped flux created by the two-sided pinching of field 
lines toward Rx- In the presence of non-zero resistivity 77, field diffusion will occur out of 
the X-region. This diffusion must be offset by fluid advection from both the exterior and 
the interior of Rx- 

In the exterior, disk inflow is induced both from back-reaction from the X-wind (if 
present) and by viscous inflow from the disk proper, with the level of kinematic viscosity v 
dictating the disk accretion-rate Mr>. If the viscosity v is only comparable to the resistivity rj, 
the disk inflow is too weak relative to diffusive penetration to produce a good inward pinch, a 
failure not conducive to the generation of X-winds. To produce sufficiently outward-bending 
field lines from viscous accretion, it is necessary to have r\ smaller than v (by a factor of 
roughly the disk-aspect ratio z /zu <^ 1; see Lubow, Papalozou, & Pringle 1994 and Shu et al. 
2007). For the specific problem of an accretion disk interacting with a stellar magnetosphere, 
the simulations of Romanova et al. (2008) demonstrate explicitly that the condition r] <C v 
is indeed the crucial ingredient to achieving an X-type magnetic configuration and thereby 
generating an X-wind. 

In the interior, a funnel flow involving closed field lines links the star to the disk. 
If the stellar rotation rate is not chosen to be the same as the angular speed ~ Qx 
of the parts of the disk to which the filed lines are rooted, rapid transients are induced. 
Some early simulations, which started with slowly rotating stellar magnetospheres linked to 
rapidily rotating disks in the initial state, managed to obtain temporary X-type magnetic 
configurations via such transients (Hayashi et al. 1996, Miller & Stone 1997, Kiiker et 
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al. 2003). The sudden removal of angular momentum from the disk by magnetic torques 
creates a dynamic onrush of material toward the star, generating a temporary magnetic 
pinch and outflow. In these studies, finite resistivity is sometimes included in the code, but 
this resistivity plays no important role because the magnetic diffusion is slow compared to 
the fast inward flow caused by the large disequilibrium of the initial state. 

Goodson et al. (1997, 1999), Goodson & Winglee (1999), and Romanova et al. (2002, 
2003, 2004a, 2005, 2008; hereafter R02, R03, R04, R05, R08) pioneered the incorporation of 
more realistically rotating stellar magnetospheres, with the later papers of the latter group 
and Goodson et al. (1997, 1999) using controlled levels of resistivity and viscosity. In 
particular, R02 and R03 performed axisymmetric and non-axisymmetric calculations that 
start with an unperturbed, aligned and tilted, stellar dipole field threading a circumstellar 
disk that begins to accrete slowly via a postulated disk viscosity v. The only resistivity 
in the problem is numerical and gives an effective value for r\ comparable to v. Although 
no X-winds arose as a consequence, these authors did show that in steady-state the radius 
R co where the magnetosphere corotates with the disk is close to the disk truncation radius 
(denoted by them as the stellar magnetopause i? m ), with the best runs having R co /R m = 
1.2-1.3 (see also Long et al. 2005). Using the same basic configuration, but rotating the 
star more quickly so that corotation is reached interior to the truncation radius - the so- 
called "propeller regime" (Illarionov & Sunyaev 1975) - R05 found strong disk outflows to 
be possible if v is several times larger than rj (see also Ustyugova et al. 2006; hereafter U06). 
The simulations of R04 discovered "magnetic towers" in the stellar corona near the rotation 
axis as first proposed by Draine (1983) and Lynden-Bell (1996). With realistic levels of 
coronal density, such towers do not carry much matter and cannot explain YSO jets (Long 
et al. 2005; see also Figs. 4 & 5 of Allen et al. 2003b for the appearance of the "magnetic 
tower" phenomenon in the context of the collapse of a rotating, magnetized, molecular-cloud 
core). Of specific interest to our present paper, the lowering of the coronal densities assumed 
by R04 by two orders of magnitude in the simulations of Long et al. (2005) made much 
easier the opening of stellar field lines, and led the way to the driving of powerful YSO jets 
in an X-type magnetic configuration by R08 when r\ <C v and R co m R m Rx- 

The main remaining differences between ideal X-wind theory and numerical simulations 
concerns whether outflows can be steady. Figure 1 shows that field lines dead to inflow or 
outflow (black) separate the funnel flow (red) from the X-wind (blue). Some of the dead- zone 
field lines are closed and link the disk to the star; some are open and "joined at infinity" 
line by line to the magnetic field contained in the X-wind. The reversal of field direction 
across the separatrix between the open dead-zone fields and the open X-wind fields of ideal 
MHD would be unstable to reconnection events in the presence of finite resistivity. These 
reconnection events, which are likely to be episodic (Aly & Kuijpers 1990), would create a 
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fluctuating X-wind (Shu et al. 1997; Romonova et al. 1998; Uzdensky, Konigl, & Litwin 
2002), and may underlie the outbursting behavior found in the simulations of Goodson et 
al. (1997, 1999), R04, R05, U06, and R08. The steady assumption of ideal X-wind theory is 
then made for analytic simplicity, and can at best represent only the time-average of outflows 
that are time-variable and quasi-periodic in reality. 

2. Generalized X-Wind Model 

To date, most of the formal developments of X-wind theory surmise that the magnetic 
field configuration of the star in the absence of any interactions with the disk is a pure dipole. 
But this is a dubious simplification because there is no reason to expect that the dynamo 
action in slowly rotating, fully convective objects (such as T Tauri stars) will create a highly 
organized pattern of surface magnetic fields. Indeed, polarization studies indicate that an 
organized dipole component cannot be dominant on the surfaces of most T Tauri stars 
(Valenti & Johns-Krull 2004). Moreover, the observed covering fraction of "hot spots" due 
to funnel flows is typically much smaller (0.1 - 1% of the stellar surface) than the predictions 
from dipole models (Johns-Krull & Gafford 2002, hereafter JG02, and references therein; see 
also Gregory et al. 2006a and Jardine et al. 2008). Recent numerical simulations by Chabrier 
& Kiiker (2006) similarly predict that in fully convective objects, the a 2 dynamo effect should 
produce large-scale non-dipolar magnetic fields, dominated by higher-order multipoles. 

Numerical simulations have kept pace with the observational developments. Thus, Ro- 
manova et al. (2004b) have performed 3-D MHD simulations of the funnel flows that result 
from inclined stellar dipoles, and Long et al. (2007, 2008) have included the further super- 
position of quadrupole contributions to such 3-D calculations. Relative to the observational 
issues raised by JG02 regarding funnel flows and their hot spots, the most important find- 
ing of these 3-D results is that the central portions of the hot spots are hotter than their 
peripheries, yielding the possibility that the total hot-spot covering fraction may have been 
under-estimated by the UV observations. 

Fortunately, the general validity of the X-wind model does not depend on detailed 
assumptions such as the dipole approximation. The important features are (a) the existence 
of trapped flux $ t in the X-region, and (b) the assumption of disk locking. Let us see how 
these two ideas lead to relationships between physical parameters of the system that agree 
better with observations of the funnel flows of T Tauri stars than other predictions (JG02). 
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2.1. Models with Trapped Flux 

We adopt cylindrical coordinates (w, ip, z) with origin at the center of the star. Suppose 
that the pinch toward the X-point traps an amount of magnetic flux equal to $i in an 
equatorial ring with radius w = Rx- In the immediate vicinity of the X-region, the trapped 
flux must be nearly force-free above and below the midplane (which is not force- free). In 
the cold limit where the X-region corresponds essentially to the single radius Rx, the field 
configuration achieves the shape of a complete fan, and a third of all the trapped field 
lines bow sufficiently outward or inward to launch an X-wind or a funnel flow, carrying 
respectively a fraction / and (1 — /) of the disk accretion rate Md- The fraction / of matter 
that physically climbs onto open field lines in the X-wind depends on how diffusive effects load 
different field lines in the X-region. Lacking such detailed knowledge, Shu et al. proposed 
guessing that / ~ 1/3 on the basis that 1/3 of the trapped field lines in the X-region bend 
sufncently outward to become opened as wind streamlines. To provide the most general 
set of expressions admissible within our analysis, in equations (1)-(13) below we assume 
that the flux fractions participating in the funnel flow and X-wind are each 1/3, as imposed 
mathematically by the cold limit of our analysis, but, in the absence of further constraints, 
keep the corresponding wind mass fraction / a free parameter (see further discussion at the 
end of §2.2). 

The field lines that are loaded with infalling gas in the funnel flow are the same ones 
that carry the gas all the way to the star. Let Fh be the fraction of the surface area 2nRl of 
the upper hemisphere covered by hot spots of mean field-strength B^. Then the magnetic 
flux connected to hot spots, Ff l (27rRl)Bf l , must be equal to the same 1/3 of the trapped flux 
$t (defined for convenience to be positive) that links to the base of the funnel flow in the 
upper surface of the X-region: 

F h (2irRl)B h = i$ t . (1) 

We will assume that the lower hemisphere looks the same as the upper hemisphere except 
with reversed directions of the magnetic field. 

In order to produce a self-consistent MHD flow, the trapped flux must be able to produce 
a time rate of change of the ^-component of the angular momentum in the funnel flow equal 
to 

(l-f)M D (l-J*)R x Q x , (2) 

where R x Qx is the angular momentum per unit mass of material orbiting the star at the 
inner disk edge, and J* is the average fraction of this specific angular momentum that lands 
on the star (a conserved quantity carried partly by matter and partly by the magnetic torques 
of the field). An analogous argument must apply to the torque supplied to the X-wind by 
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the flux $t/3 contained in it; namely, the wind torque must supply a time-rate of change of 
the angular-momentum equal to 

fM D (J w -l)R 2 x Q x . (3) 

Thus, angular momentum enters the X-region at a rate MdR x Qx through disk accretion, 
while angular momentum leaves the X-region at a rate f Me,J w R 2 x VLx through the outflow in 
the wind, and at a rate (1 — f)Mr> J*R x Qx by accretion through the funnel flow. In addition, 
there is a viscous toque T > exerted by disk matter interior to the X-region on disk matter 
exterior to it. In steady state, the single term representing angular momentum entering the 
X-region must be balanced by the sum of the three terms leaving it. If we solve the resulting 
equation for the fraction / of the disk accretion rate that is carried in the X-wind, we get 
(see Shu et al. 1994a) 

/ = ^^, (4) 

where r = T/M D R X Q X > is the dimensionless viscous torque acting across a circle beyond 
Rx- 

Observationally, it is possible to measure J w in principle from the mean terminal velocity 
v w reached at infinity by the material in the X-wind. For a cold flow, conservation of energy 
in the corotating frame (Jacobi's constant) implies that (see Shu et al. 1994a): 



v w = Rx^xV^Jw - 3. (5) 

With the guesses that J* ~ for a small, slowly rotating star, and r ps 0, equations (jl]) and 
dSJ) yield J w ps 1/ f ps 3 and v w ps y/3Rx^x, which are in rough agreement with observations 
of YSO winds and jets (Shang et al. 2002, 2004; Pyo et al. 2006). To obtain the standoff 
distance Rx for the stellar magnetopause, we need a model that gives the trapped flux $ 4 
as a function of distance Rx from the center of the stellar multipole distribution (see §§2.2 
and 4.1). 



2.2. Connection with Dynamics of X-Wind 

To obtain the X-point location in the meridional plane, Rx, it is necessary (as in OS95) 
to input the dynamics of the X-wind because, by assumption in the simplest model, it is the 
X-wind back torque that is trying to drive the inner edge of the disk inward. The funnel flow 
field lines, being strong, then provide whatever (small) azimuthal component B v is needed 
to yield a back torque that holds the inner edge in equilibrium at Rx- This assumption 
ignores any role for a magnetic tower that might torque down the spin of the star (e.g., R04, 
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Long et al. 2005), and it is an open question in our opinion whether real YSO coronae can 
be sufficiently dense to make such braking important. We assume that 1/3 of the trapped 
flux $j drives an X-wind given by (see eq. 2.3 of OS95 and 3.10b of Shu et al. 1994b): 

U t = 2-KBf 1/2 {GM^MlR\) l /\ (6) 

where (3 is a dimensionless (inverse mass-loading) parameter that measures the ratio of 
magnetic field to mass flux in the a frame that corotates with the footpoint of a field line 
labeled by the streamline ijj averaged over X-wind streamlines, if) = to tp = 1. That is, 
B = f3(ip)p\i with /3 = J /3(ip) dip. The normalization for the streamfunction is the mass-loss 
rate in the wind, M w , which is a fraction / of the disk accretion rate Mp. The quantity 
{GM^M D R X ) 1 ^ in equation dSJ) then provides the correct units when GM», Rx, and Mb 
are taken to be the fundamental dimensional quantities of the physical problem. 

Setting <E>t = 6nBhFhRl from equation (1), we obtain from equation (jSJ): 



D \3/4 

F h B h = Pf' 2 B m ( ^ , (7) 



where we have defined a fiducial field strength, 



GMJvIf 
Rl 




We emphasize that the trapped-flux fractions, 1/3 for the funnel flow, 1/3 for the dead 
zone, and 1/3 for the X-wind, result from assuming that the thermal speed near the inner disk 
edge is small in comparison to its Kepler speed. This assumption allows us to shrink the X- 
region to a mathematical point in the meridional plane, with the disk having infinitiesimal 
vertical thickness in the same approximation. How the vanishingly small diffusivities v 
and r] then load field lines is not an addressable question in this limiting procedure, so 
we have to make a guess for this process. For the somewhat ad-hoc loading law in the 
X-wind, (3 = (3(3/3/2) (1 - ipy 1/3 , Table 3 of Cai et al. (2008) yields the identifications: 
J w = 2.64, 4.36, and 6.20, respectively, for /3 — 1, 2, and 3. With the further assumption that 
f — 1/J W — 1/3 made for simplicity, we can then interpolate to obtain the identification 
that (3 = 1.21. 
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2.3. Models with Disk Locking 

We now impose the disk locking condition, 

With the stellar rotation rate given by the above, we obtain from equations ([1]) and ([6]) the 
desired general relationship: 

F h RlB h = $ t = pf/^GM^Mn/tt*) 1 / 2 . (10) 

Equation ffTUj) . the exact analogue of a scaling relationship first derived by JG02, encapsulates 
the idea of flux trapping (see §2.5) because it relates the amount of measured flux in hot 
spots on the left-hand side to independently observable quantities of the system on the right- 
hand side, without any assumptions about the multipolar character of the stellar magnetic 
field ultimately responsible for the funnel flow. 



2.4. Dipole Models 

For the case of an undisturbed magnetic configuration for the star with a pure dipole 
of moment //*, OS95 found the trapped stellar flux required for steady state to be 

*, = ^. (id 

tlx 

In the limit of vanishingly small resistivity, we can picture the trapped flux in equation (1111) 
to originate by the disk accretion flow sweeping all the exterior dipole flux from oo to Rx, 
2ti[i*/Rx, into the X- region, plus another half that amount, tt^/Rx, resulting from the 
interior dipole flux from 2Rx/3 to Rx being pushed into the X- region by the back reaction 
to the outward transport of angular momentum in the funnel flow. 

Putting the expression (iTTj) into equation (jSJ), we may solve for Rx to obtain 

^•yy" <i2) 

where at is given by 

at= {2W^) ' (13) 
an expression first derived by OS95. From the early computations of Najita & Shu (1994), 
OS95 estimated (3 = 1 would be needed to drive an X-wind with / = 1/3, giving a t = 
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0.921. (OS95 actually obtained 0.923 because eq. [TT] represents a slight rounding of their 
numerical result for dipole funnel flows.) With the more refined estimate (3 = 1.21 for the 
case / = 1/3 from Cai et al. (2007) as discussed earlier, we get a t = 0.825 for the dipole 
model. More general multipole models of funnel flows do not satisfy a simple relationship 
for the magnetopause standoff distance such as equation (Tl2l) ; the dependence on the stellar 
dipole moment /i* is replaced by a more complex dependence on the field geometry on the 
stellar surface at radius R*. This is encapsulated in the discussion of §§3.3 and 4.2, where 
the surface field geometry helps to determine the quantities F h and B h . 

Although equation ffl2l is reminiscent of similar expressions in Ghosh & Lamb (1978, 
1979a, b), we emphasize that the resemblance arises mostly as a result of dimensional analysis. 
Given the three dimensional quantities GM*, /z* and Md-, there is only one combination, 
{l^i/GM^Mp) 1 / 7 , apart from a numerical coefficient a t of order unity, that yields a quantity 
Rx with the dimensions of length. The physical reasoning underlying Ghosh & Lamb's 
derivation for the disk truncation radius is quite different from those given above. Specifically, 
when steady state holds in a problem where the resistivity is small, the disk is nearly in 
mechanical equilibrium with stellar gravity balancing the centrifugal force of rotation. Thus, 
there is no unbalanced "ram pressure" that the nearly corotating magnetic field at the disk's 
inner edge needs to offset. What the stellar magnetosphere needs to do to truncate the 
disk is to transfer enough positive angular momentum outwards so as to offset the negative 
torques that are trying to drive the disk inwards in the equatorial plane. 

The important point that disk truncation in CTTSs arises, not from ram pressure effects, 
but from considerations of angular momentum transport was first made by Cameron & 
Campbell (1993). Their analysis is closest in spirit to OS95's and ours. The similarity is 
nearest in the prediction that disks are truncated at inner edges which lie close to the radius 
where the disk corotates with the central star, a phenomenon that has come to be called 
"disk locking." Both analyses also ignore the presence of any intrinsic magnetization of 
the disk. The main differences are that we suppose (1) that the disk resistivity is much 
smaller than the disk viscosity (by a factor proportional to the disk-aspect ratio; see Shu 
et al. 2007), and (2) that the viscous torque in the disk is negligible relative to the funnel- 
flow and X-wind torques, so that neither of the values of the coefficients of viscosity v and 
resistivity rj enter the final formulae in the limit of a vanishing ratio for the sound speed 
to Keplerian velocity in the X-region, except for the dependence implicit in our assumption 
for the wind-loading fraction /. In contrast, Cameron & Campbell (1993) make detailed 
assumptions concerning v and 77, which they assume are equal in turbulent circumstances, 
and they ignore the X-wind torque, equating the funnel-flow torque to (the negative of) 
the viscous torque. Moreover, instead of computing the funnel-flow torque in the funnel, 
they estimate it from the slip present between the mid-plane of the disk, assumed to rotate 
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at Keplerian speeds, and its surface, assumed to corotate at the angular speed of the star, 
for regions that are magnetically connected to the star. They parameterize this slip by an 
unknown factor 7 C 1 that characterizes the resulting azimuthal magnetic field giving rise 
to a magnetic coupleQ If 7 ~ 1, then their R m (= Rx in our notation) can be larger than 
our equation ({121) by as much as a factor of 2 in typical circumstances (see §4.3). Given 
the modeling uncertainties and the effects ignored in each exercise to obtain a tractable 
analysis, we do not consider this difference as numerically significant at the present state of 
observational tests. 

In any case, the substitution of equation ffT2l into equation ([9]) and /x* oc BhR^ results 
in a proportionality relationship specific to the dipole case, 



This is to be compared to the general expression for arbitrary stellar fields, equation ( fTOl) . 



JG02 have compared observations to the dipole model predictions by various groups, 
including that of OS95 embodied in equation (|T4|) . as well as to the multipole prediction en- 
capsulated in equation (JTUJ). To make progress, they assume that all T Tauri stars have the 
same value of Bh/fif 1 / 2 , allowing both equations (ITUj) and (TBI) to be treated as proportion- 
alities. In the restricted version of X-winds presented here, where (3 = 1.21 and / = 1/3, this 
assumption amounts to postulating that Bh has the same value in all hot spots, a supposition 
that seems to hold approximately in many observed cases. They find that the observations 
bear little correlation with all dipole predictions, including that of OS95. The scatter is 
greatly reduced and a significant correlation is found, however, when the data are compared 



In the notation of this paper, equation (3) of Cameron & Campbell (1993), which defines 7, reads 



where is the mean toroidal magnetic field at the surface of a disk of half- height Zq, and f2» and SIk 
are, respectively, the rotation rates of the star and the disk (which is Keplerian) in the neighborhood of the 
truncation radius w. To avoid violent Parker instabilities, Cameron & Parker require \B^\ to be not much 
larger than B z , whereas their analysis suggests that VLk and SI, can differ by amounts that make w\Qk — 
considerably larger than the thermal speed a in the midplane of the disk. They consider two cases, where 
rj ~ azg/vu and r\ ~ azQ (cf. the unnumbered eq. after eq. 3 and eq. 7 of their paper); in both cases, their 
equation (3) requires 7 <C 1 for disks that are spatially thin, i.e., for which zo/w <C 1. Only a stretch of all 
parameter choices could make 7 ~ 1. 



B h Rl oc {GM^M^ 2 ^ 6 . 



(14) 



2.5. Analysis by JG02 
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to the relationship of equation fflOl) . where the dipole approximation is dropped. Thus, the 
generalized X-wind model for accretion flows is indeed compatible with current observations. 

We note that the mean hot spot field strength, B^, is also in principle measurable by 
Zeeman effects. Thus, if Pf 1 ^ 2 is a theoretically computable coefficient, then not only can the 
slope of a log-log plot of the variables on the left- and right-hand sides of equation (fit)]) be 
checked, but so can its intercept. Indeed, Bh and Fh have now been simultaneously measured 
in at least two T Tauri stars (V2129 Oph, Donati et al. 2007; BP Tau, Donati et al. 2008); 
we compare the data to our theory in §4.3 and find reasonable agreement. 

3. Numerical Calculations 

We wish to pursue the consequences of the above theoretical discussion. In particular, we 
wish to test whether flux-trapping at the X-point can lead to accretion funnel flows even when 
the unperturbed stellar field is not a pure dipole. Our technique is to use direct numerical 
calculation, analogous to that carried out by OS95 for the dipolar case. Specifically, we will 
model time-independent, axisymmetric accretion flows where the unperturbed field on the 
stellar surface is composed of higher-order multipoles. We state at the outset that we adopt 
the simplification of OS95 that the sub-Alfvenic conditions of the funnel flow and the dead 
zones surrounding it can be approximated with vacuum fields plus current sheets, rather 
than treating the problem by the more accurate, but also considerably more complex, Grad- 
Shafranov equation (see, e.g., Cai et al. 2008). Even with this simplification and the further 
assumption of axial symmetry, the removal of the dipole constraint means that our particular 
choices of the stellar field configuration will be neither unique nor universally applicable. We 
expect the actual field to differ from star to star, and to vary with time on any given star. 
The goal of our study is more limited: to demonstrate explicitly how it is possible to satisfy 
the two primary observational constraints that the net surface field polarization in T Tauri 
stars is small, and that the accretion hot spots cover only a small part of the stellar surface. 

3.1. Multipole Model 

The dipole model calculations for the inner (funnel-flow) solution have been discussed 
in detail by OS95. Here we describe only our new multipole calculations for the funnel 
flow, while drawing attention to the salient points of departure from the dipole model. The 
particular stellar multipole field we choose is discussed in the next section. 

We assume an unperturbed stellar field that is some multipole aligned with the stel- 
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lar rotation axis. We also assume that the poloidal component of the final field (i.e., the 
perturbed field after accounting for interactions with the disk and flux-trapping) is in a 
vacuum- field potential configuration (i.e., is curl-free), including the field lines loaded with 
accreting matter. This condition is approximately satisfied provided (i) the Alven Mach 
number of the flow, M.A-, satisfies M.a *C 1 everywhere, and {%%) the net fractional specific 
angular momentum carried to the star in the funnel flow, X, satisfies | J* | <C 1. The first 
condition means that the magnetic field energy dominates the matter kinetic energy even 
on the field lines participating in the funnel flow. This is certainly true at both the X-point 
pinch and near the surface of the star, given the high field strengths in those regions. With 
the condition satisfied at both ends, we (like OS95) plausibly assume it is also true every- 
where between. The second condition means that any stellar spin-up or spin-down, due to 
the specific angular momentum carried on to the star by the accretion flow, is very slow - of 
order the accretion timescale - as is observed (see discussion of equation (2) in §2). These 
considerations justify our (and OS95's) assumption of a curl-free poloidal field, V x B p = 0; 
i.e., the associated magnetic flux <D satisfies the vacuum-field condition, 



OS95 adopt the agnostic choice J* = to keep a small star slowly rotating. In this 
case, with M. A <c 1, the field lines in the funnel flow must have a small but non-zero 
azimuthal component. Combined with the poloidal field, this produces field lines that are 
trailing spirals, leading to an outward transfer of angular momentum, extracted from the 
inflowing material in the funnel flow by the resulting magnetic torques, and transferred 
back to the footpoints of the funnel flow field in the X-region of the disk. Because the 
azimuthal component is small compared to the poloidal field in the funnel flow, the overall 
field geometry interior to the X-point is well represented by the poloidal field B p alone. 
The latter is calculated by solving equation (fl5|) subject to appropriate boundary conditions 
(described below). The small azimuthal field B v can be found a posteriori, through the 
details of how the field lines are loaded with matter. 

Calculating the axisymmetric B p then reduces to a 2-D problem in the meridional plane. 
Since we wish to solve only the inner problem (funnel-flow accretion), our calculations are 
confined to the region enclosed by the boundaries z=0 (axis along the equatorial plane, 
i.e., along the plane of the disk), w=0 (rotation axis of the star and symmetry axis of 
the stellar multipole, perpendicular to the disk plane), z=z max (upper boundary, parallel 
to the equatorial plane, of our modeled region) and Z — Z w (ijJ ) (locus of the uppermost, i.e., 
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innermost, wind streamline) (see Fig. 1). The corresponding boundary conditions are: 



$ -> $ m as [ro 2 + z 2 ] -> (16a) 

$ -> $ f as [(uj - i? x ) 2 + z 2 ] -> (16b) 

$ = $„, = constant on z = z w (vj) (16c) 

$ = $. a ,(ro/a7 max ) 2 onz = z max (16d) 

$ = $ t on 2 = for ^ < ro < R x (16e) 

d$/<9z = on ^ = for < w < R k (16f) 

$ = on t*7 = 0. (16g) 



In the above, $> m (zj,z) is the unperturbed multipole flux; $ f = $ t #/7T is the flux 
function for a fan of poloidal field over the angles < $ < tc measured in the meridional 
plane with vertex at the X-point: *& = arctan [z/(w — Rx)]- The quantity $ t is the total 
amount of flux trapped between •& = and 7r; $> w = $t/3 is the flux adjacent to the 
uppermost wind streamline. The radius Rk is the position of the "kink point" where stellar 
magnetic field lines make a transition from skimming horizontally over the equatorial plane 
for w > Rk to threading it vertically for w < R k , because the gas falls (in a trailing spiral) 
from the magnetic fan emanating from the X-point toward the star along the lowermost 
funnel-flow field lines that have been pulled into the inner edge of the accretion disk in 
back reaction to the torques exerted by the funnel flow. A more detailed explanation of the 
boundary conditions (16a)-(16b) is supplied in the Appendix, together with a prescription of 
how they are used in conjunction with the governing PDE (ITS]) to find the magnetospheric 
configuration interior to the X-wind (see Fig. 1). 



3.2. Specification of Stellar Multipole Field 

Written in terms of a scalar potential a poloidal field satisfies (Chandrasekhar 1961): 



B = V x 



x r 



v r 

Imposing the curl- free condition V x B p =0 reduces this to Laplace's equation for ty/r: 

v^)=o. 

r / 



(17) 
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In a multipole expansion of spherical harmonics, this has the solution: 

oo oo 

¥(r, n,<p) = E ^ r ^ ^™ r ' +1 ^' ( 19 ) 

1=0 m=0 

where /i = cos 8. We look for physical solutions wherein the strength of each multipole 
component of the unperturbed field decreases with increasing radius (distance from the 
star). We also wish to consider only azimuthally symmetric solutions independent of if, i.e., 
with azimuthal quantum number m—0. With no 92-dependence, the associated flux is 

* m (r,„) = (20 ) 

Putting the above equations together, the unperturbed multipole flux can be expressed as 

$ m (r,^) = -(l- / . 2 )f;^r-'^M, (21) 

1=0 ^ 

where Pi(n) are the usual Legendre polynomials. The function $ m is known globally once it 
is fixed at any given radius r = R because the coefficients in equation (121]) are 

A _ (21 + 1) Rl f +1 i [Q m {R,n) 



J [J i*', df] Pi(ji) dp. (22) 



Jardine, Cameron, & Donati (2002) use a similar method to represent the surface magnetic 
fields of real stars derived from Zeeman tomography. 

We elect to define the multipole flux $ m at the stellar surface, r = i?* = 1 for the 
present normalization. To choose the surface flux, we note that the radial magnetic field on 
the surface is given by 

We impose the conditions that the radial surface field is an odd function of /i: B r (R*, /i) = 
—B r (R*, — /i), so that it vanishes at the equator: B r (R*,fi) = for \i = 0. By equation f l23l) . 
the same conditions must hold for $^ = <9$ m (-R*, fi)/dfi. Up to an arbitrary multiplicative 
constant, a generic odd function that satisfies these conditions is 

^(l^)=/i 2p+1 (l-^ 2 ) 9 , (24) 

which reaches maximum/minimum (corresponding to maximum inward/outward directed 
B r ) at the latitudes /io = ±a/ (2p + l)/(2p + 2q + 1). This form allows one to control both 
the compactness and the ingress/egress latitudes of the surface radial field, and hence of the 



-17- 



final hot spots. Multiple loops can be created by a superposition of such states with different 
exponents p and q, and/or multiplying by a suitable oscillating even function in fi. 

For demonstration purposes, in this paper we use one such surface field: 

$; n (l,/i) = cos(80)//(l -ii 2 ). (25) 
Expanding cos(8#) in terms of /i = cos 6, we obtain 

$' m (l, 0) = [1 " 32^ 2 (1 - /x 2 ) (1 - 2^ 2 ) 2 ] /i (1 - /x 2 ) (26) 

With $^(l,/i) specified, and the normalized radial scale R* = 1, we have $ m (l,/i) 
= f $^(1, /i) d/i. (The constant of integration can be set to zero without affecting the B- 
field). This expression is inserted into equation (22) to solve for the coefficients Ai, which 
are then employed in equation (21) to calculate $ m (r, fi) globally. Note that equations 
fT24|) - fr26|) express the non-dimensional flux on the stellar surface with an arbitrary scaling. 
The dimensional value of the corresponding surface magnetic field depends on the stellar 
parameters, the accretion rate to be supported, and the details of how matter is loaded onto 
the field lines. This physical scaling was discussed in §2.2 (see equations (7) and (8)), and 
is made concrete for our specific models in §4.2. 

Our choice of $ m , obtained from equation (1231) . corresponds to 6 multipole components: 
I = 1, 3, 5, 7, 9, 11, where I = 1 is the dipole part. The corresponding Ai, relative to the 
dipole component A\, are given in Table 1. Figure 2 shows the resulting B r and Bg on the 
stellar surface and Bg as a function of radius in the equatorial plane (note that B r = in 
the equatorial plane). Figure 3 shows the unperturbed magnetic field lines over the entire 
computational domain. The salient features of our chosen field are as follows. 

The Ai are all roughly equal within a factor of a few, implying, via equation ( T2T1) with 
r = i2* = 1, that the various components of the multipole flux are all roughly comparable 
(within an order of magnitude) on the stellar surface (the corresponding radial magnetic 
field on the surface, B r (r=l) is dominated by the higher orders, with a very weak dipole 
component). The corresponding surface field geometry is highly loopy, as implied by the 
oscillatory shape of B r and Bg in Figure 2 and explicitly illustrated in Figure 3. Thus the 
net polarization of the starlight derived from the surface field will be low. Moreover, since 
the radial field is confined to relatively narrow latitudinal bands, the self-consistent location 
of the X-point must be found so as to capture a significant fraction of the magnetic flux 
in these latitudinal bands; otherwise, there will not be enough trapped flux. Thus, the 
hot-spot covering fraction will also be low. Hence, highly loopy field configurations are 
automatically consistent with the observed surface properties of CTTSs. 
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Moving away from the star, the higher order multipoles rapidly weaken, with Bi(r) 
— > r^^ +2 ^; the loopiness of the field decreases; and the purely dipole component finally 
dominates at large distances (at > 5R* for our adopted field; see Fig. 3). We can therefore 
expect the funnel flow solution near the X-point (i.e., the values of $ t /$ mx and Rk/Rx) to 
resemble the OS95 dipole solution when Rx/R* is sufficiently large. However, we emphasize 
that, independent of the value of Rx/R*, the important consideration is the trapped flux, 
and not the detailed multipole character of the stellar field. 

It is noteworthy that the unperturbed field varies non-monotonically in both strength 
and direction with radius, unlike in the dipole case. While each of the multipole components 
declines monotonically with radius, the vectorial directions of all the components are not 
the same. Consequently, the multipole fields cancel out more in some regions than in others, 
yielding a summed field that oscillates in strength and direction with radius and is null at 
some radii, specifically, at 1.015/2* and 3.25/?* (Figs. 2, 3). The summed field <3> m (r, jj) 
finally decreases smoothly to zero at large distances, once the dipole dominates. For the 
pure dipole case, OS95 normalized by Rx to find a unique solution (in terms of $t/$dx 
and Rk/Rx) independent of Rx/R*- This strategy is possible since the unperturbed field 
looks the same geometrically at all radii for a dipole (or for any other single higher order 
multipole), and the location of the stellar surface is immaterial for the scaled solution. In our 
case, the location of the stellar surface is fixed by our choice of the normalization length for 
the specification of different multipoles. As a consequence, the funnel flow solution (& t /&mx 
and Rk/Rx) for a specified surface field configuration composed of a mixture of mutipoles 
cannot adapt unchanged to a stellar radius R* that is chosen arbitrarily relative to Rx- 

However, once a funnel flow solution is found for a given Rx, the stellar surface can be 
redrawn at an arbitrary location to represent the solution for a different surface field than the 
one originally specified. This is possible because, while the initial unperturbed field is fixed 
globally by its specified configuration on the stellar surface, our actual numerical calculation 
of the solution from this depends only on the location of Rx relative to the origin and not 
with respect to the stellar surface (i.e., the stellar surface is not used as a boundary condition, 
the origin is). Thus, while the derived solution corresponds to the originally specified surface 
field if the surface is drawn at R*=l, it also corresponds to a different surface configuration 
for R*^l (with the relative surface field strength contributions of the n multipoles in the new 
configuration being R~ 3 : R~ 5 : ... : /?™~ 3 times their originally specified values). Obviously, 
the solution for a given Rx will correspond to a different Rx/R* for R*^l compared to 
R*=l. However, since the solution fixes the amount of trapped flux, the quantity F h B h R\ 
- the funnel flow magnetic flux on the stellar surface - remains constant independent of the 
R* location (see equation [1]). In dimensional units, if R* is held constant (and M* and 
Md are fixed, so B norm remains unchanged), then the quantity F^Bh is uniquely determined 
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by Rx/R* (see equation [7]), with each new redrawing of the stellar surface equivalent to a 
different Rx/R*- What changes, given a solution at fixed Rx, is the spot covering fraction 
F h and corresponding B h (with F h increasing as Rx / R* decreases, since the funnel broadens 
with distance from the star, and B h thus decreasing). We shall use such redrawings to 
compare the dimensional results of our calculations to observations, in §4.2 and 4.3. 



3.3. Covering Fraction of Hot Spot 

For axisymmetric configurations, the covering fraction Fh of the hot spot in the upper 
hemisphere is given by 

r-e-2 

F h = / sin0d0 = ^ -/i 2 , (27) 



where /i = cos 9 and 9i and 6* 2 are, respectively, the upper and lower colatitudes of the 
hot-spot (annulus). The numerical values of fii and /z 2 are given by the flux tubes that link 
them to the upper and lower surfaces of the trapped flux in the X-region, i.e., 

$(r = R*, f i 1 ) = $(r = R x ,# = ir); $(r = R*, /x 2 ) = $(r = Rx, $ = 2vr/3). (28) 

When Rx 3> -R*, the hot spot is generally small, and we may approximate the flux difference 
on the stellar surface as arising from the unperturbed stellar field, $(.??*, /z 2 ) — $(i?*,/ii) ~ 
$' m (i?*, /x/ l )(/ii — /i 2 ), where fih = (a*i + A^)/^ gives the mean latitudinal location of the hot 
spot, and Q' m (R*,fi) in unnormalized dimensionless form is the function fT2B|) . Moreover, the 
difference in $ at Rx between the angles d = 27r/3 and 7r is simply 1/3 of the total trapped 
flux $ t . Thus, when the covering fraction is small, we have the approximation, 

F h = f i 1 - f i 2 ^—^ -, (29) 

with it being unimportant in the ratio on the right-hand side whether we use scaled dimen- 
sional values for $ t and &' m (l, fJ-h) or unsealed nondimensional values. 



4. Numerical Results 

4.1. Steady State Solutions 

We present dimensionless solutions for the X-point located at four representative radii: 
Rx = 10-R*, 7.5i?*, and 2R*. In other words, we assume here that the dimensional 
parameters of the problem are such that the equilibrium value Rx has the chosen ratio 
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relative to i?*, before we consider in §4.2 the requisite value of the mean hot-spot field 
strength needed to accomplish this feat. The results are as follows. 

R x = 10 R* '■ The final field configuration for this case is plotted in Figure 4. The require- 
ments are well met for a solution in which we may expect a sub-Alfvenic funnel flow to 
occur in the manner described by OS 95: the uppermost funnel flow field line just skims 
the equipotential curve near the X-point and never rises above this curve, ensuring steady 
accretion onto the star, and the lowest funnel flow field line subtends an angle of 120° to 
the disk as it exits at the kink point. The accreting material lands on the star at a mean 
colatitude of 6 = 42.0°. The solution corresponds to $t/$ mx = 1.55 and Rk/Rx = 0.70. 
This is almost exactly the solution OS95 found for the dipole case: $ t /$d x = 1.5 and Rk/Rx 
= 0.74. As discussed above, this is to be expected, since the dipole component of the field 
vastly dominates at 10-R* even in our multipole model. The situation is dramatically dif- 
ferent close to the star, however. Comparing OS95's dipole solution plotted in Figure 1 to 
our multipole one in Figure 4, we see that the flow in our model is squeezed into a much 
narrower funnel close to the star, as the higher order multipoles begin to dominate. Long 
et al. (2008) noticed a similar trend in their 3-D numerical simulations when the stellar 
magnetic configuration is complex. In fact, the resolution of the inkjet printer used to make 
this figure is too low to show the true smallness of the hot spot; our hot spot covering frac- 
tion as computed by equation ff29l is only 0.065% (latitudinal width of only AO = O.O56°)0, 



compared to the much larger ~ 3% for a pure dipole with Rx = 10R* (see OS95)13 

Rx = 7.5-R* : The converged solution for this case is plotted in Fig. 5. It corresponds 
to $t/$ mx = 1.71 and Rk/Rx = 0.61. The flow lands on the stellar surface at a mean 
colatitude of = 42.1°, with a covering fraction of 0.093% (AO = 0.08°), compared to ~ 4% 
for a pure dipole with R x = 7.5-R*. The solution parameters are close to the 10.R* case, but 
not identical, for the following reason. While the dipole component utterly dominates at 
7.5.R*, the quadrupole starts to be non-negligible at smaller radii, with a polarity opposite 
that of the dipole. This is evident in Fig. 5, where we see that Rk, the inner edge of the 



2 In this and all other solutions below, the covering fraction computed from equation (29) is checked 
against that directly implied by the position of the first and last funnel flow field line on the stellar surface in 
the numerical models plotted in Figs. 4-7. They agree in to within 0.01°; the analytic values from equation 
(29) are quoted as the more accurate ones, given the finite spatial resolution of the numerical calculations. 

3 Since our calculations assume axial symmetry, strictly speaking, we should speak in terms of "hot rings" 
rather than "hot spots." However, the latter is the common usage among observers in the field, and we will 
follow this usage under the assumption that the qualitative concepts of the problem important to this paper 
will not change - as is indeed demonstrated, for example, by the simulations of Long et al. (2007) - if we 
were to perform a calculation that dropped the assumption of axial symmetry. 
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inward excursion of the disk, is not very far from the magnetic null point at 3.25/?*, where 
the dipole field is completely cancelled by higher orders. Thus the net flux swept into the 
X-region increases slower with decreasing radius than in the purely dipolar case, requiring 
a larger inward extension of the disk, i.e., smaller R k /R x , to trap sufficient flux than for a 
pure dipole. 

The above behavior raises the interesting question: What happens as Rx continues to 
approach the null point at 3.25/?* from radii exterior to it? In the 7.5/?* case above, enough 
flux resides between the null point and Rx to provide a good solution, despite the increasing 
importance of the quadrupolar component of opposite polarity to the dipole one. At some 
stage, however, this will no longer be true, and flux from radii interior to the null point will 
start getting swept into the X-region exterior to it. Since the field reverses at the null point, 
the net effect will be the reduction of trapped flux in the X-region. A converged solution 
will then arise only when enough flux from radii interior to the null point has been swept 
in to offset this cancellation and once again provide sufficient trapped flux. The field lines 
supporting the funnel flow must consist now entirely of ones that originate at low latitudes, 
below the second loop, instead of above it as in the 7.5/?* case (see Fig. 3). The result is 
a latitudinal change of 30° in the hot-spot location (and a reversal in the polarity of the 
radial magnetic field in the hot spot compared to the previous location). In other words, at 
some critical juncture, the hot-spot position must jump significantly for a small shift in Rx, 
or equivalently, for a small change in the accretion rate. This is illustrated in the following 
solution. 

R x = 5/?* : The final field configuration for this case is plotted in Fig. 6. The solution 
parameters are <3>t/<I> m x = -25.0 and Rk/Rx = 0.247. The flow lands at a mean colatitude of 
9 = 72.4° on the stellar surface, with a covering fraction of 2.74% (A6 = 1.65°), compared 
to ~ 6% for a pure dipole with R x = 5/?*. These numbers validate our thought experiment 
above. The cancellation of flux in the X-region by field of opposite polarities requires far more 
flux from radii interior to the null point to be swept in to finally acquire enough trapped flux; 
consequently the inward excursion of the disk is much larger (Rk/Rx much smaller), and 
the ratio of trapped to unperturbed flux much higher, than in the 7.5/?* and 10/?* 

cases above. The negative sign of <3> t /<3> mx denotes that the net trapped flux is oppositely 
directed to the original unperturbed flux at the X-point. As a result, we see that the hot 
spot has now jumped to below the second loop, as predicted. 

The latter result resembles the conclusion reached observationally by Gregory et al. 
(2006) who deduced funnel streams flowing close to the equatorial plane when the central 
star has a complex magnetic topology. The 3-D simulations of Long et al. (2007, 2008) also 
obtained matter to accrete onto hot spots close to the equatorial plane in such circumstances. 
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Rapid changes in hot spot location are indeed observed on CTTSs, a phenomenon 
impossible with a purely dipolar field (or any other single higher order multipole). The 
discussion and calculation given above illustrate that such changes can only be explained by 
stellar fields that are complicated structures consisting of a mixture of multipoles. Non-axial 
symmetry and time-dependence would add extra richness to this phenomenon. 

Rx = 2R* '■ The converged solution for this case is plotted in Figure 7. The accreting gas 
now lands on the star at a mean colatitude of 9 = 69.5° with a covering fraction of 0.30% 
(A9 = 0.18°) compared to the very large value of ~ 21% for a pure dipole with R x = 2R*. 
The solution parameters are now <3>i/<I> mx = 40.0 and Rk/Rx = 0.87; as in the 5i?* case, 
these are significantly different from the solutions at 7.5 and 10R*. The reason again lies 
in the oscillatory behaviour of the magnetic field. The X-point traps all the flux exterior 
to Rx, as well as some fraction of the flux interior to it. Since the null point at 3.25-R* is 
exterior to Rx = 2R*, many of the field lines swept into Rx, from radii both interior and 
exterior to it, cancel out. Thus far more flux needs to be swept into the X-region (relative to 
the unperturbed flux residing there) from radii interior to it, to offset this cancellation and 
support a stable funnel flow, than at say Rx = 10-R*, where no such field turnover exists 
(either exterior to it or interior to it for radii >Rk)- This explains the much larger <3> 4 /$ mx 
in the 2R* case. At the same time, the field strength (Bg) in the equatorial plane rises much 
more rapidly, i.e., flux lines are much more tightly packed, with decreasing radius at 2/2* 
than at 10-R* (Figs. 2, 3). Thus a relatively small excursion of the disk inward of Rx = 2i?* 
is required to capture the requisite flux for a steady-state solution. This accounts for the 
smaller Rk/Rx at 2R*, in spite of the flux cancellation, than at 10-R*. The rapid increase in 
stellar field inward from Rx = 2R* implies that it would be much harder to crush the stellar 
magnetosphere entirely (i.e., bring Rx inside R* by increases of Mr,) than naive estimates 
with pure dipole formulae might suggest. 



4.2. Dimensional Field Strengths 

Typical numerical values for CTTSs are M* = 0.5 M e , R* = 2R Q , and Mo = lO~ 8 M yr _1 . 
Our range in simulated R x , 2-10.R*, then corresponds to disk-locked stellar rotation peri- 
ods of 1.3 to 14 days, consistent with the observed span in CTTS periods. These fiducial 
parameters also yield, via equation flS]), -B n orm = 26.68 G. With ft = 1.21, / = 1/3, and F h 
= 3.0xl0~ 3 , 2.74xl0- 2 , 9.3xl0" 4 and 6.5xl0~ 4 computed for the cases R x /R* = 2, 5, 7.5, 
and 10, equation ([7]) then implies that the required mean strength for the hot spot field is 
Bh — 11, 2.3, 91, and 161 kG, respectively. Because we have fixed the surface geometry, the 
absolute value of the field strength (without regard to sign), averaged over the entire stellar 
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surface, is ~60% of the hot spot field in each case. 

The measured mean absolute field strengths on CTTSs are a few kG (e.g., Johns- Krull 
2007). Formally, this suggests that R x ~ 5.R* is the expected standoff distance of the stellar 
magnetopause for the choice of field geometry and stellar and disk parameters (M* = 0.5M Q , 
R* = 2R Q , M d = 1 x 1O~ 8 M ) made here. The ratio Rx/R* ~ 5 compares well as an 
average with the measured inner edges of CTTSs by infrared spectroscopy (e.g., Carr 2008). 
However, the specifics of the field reversals in the 5.R* model gives Figure 6 an especially 
large spatial extent from R& to Rx, suggesting, perhaps, especially large X-rays flares when 
such regions undergo magnetic reconnection. It would be interesting in this regard to map 
the magnetic field configurations on the stellar surface and their loops above the surface 
relative to the locations of the hot-spots of known YSO sources of large X-ray flares. A 
pre-flare configuration that resembles Figure 6 and a post-flare configuration that resembles 
Figure 7 would be most suggestive of some of the mechanisms proposed in this paper for the 
reasons behind hot-spot jumps. 

In contrast, the B h we compute for R x = 7.5 and 10.R* are much too high, by almost 2 
orders of magnitude, to represent realistic mean absolute field strengths on CTTSs. Concur- 
rently, the hot spots in these two cases, while smaller than for a dipole, as desired, are too 
small compared to observations. Field configurations with large amounts of power in high-/ 
multipole components are unlikely to allow large magnetopause radii Rx compared to i?*, 
for the following reason. By construction, the X-region traps all flux from the outer edge of 
the disk to R^. Over this range of radii, the dipole component dominates in the equatorial 
plane for 10 and 7.5R* (with the quadrupole component making a small contribution for 
7.5i?*). On the stellar surface for these two cases, however, higher order multipoles (specif- 
ically, the 1=5 and 9 components) characterize the hot spots. In other words, the trapped 
flux that carries the funnel flow is predominantly dipolar in the X-region but overwhelmingly 
a superposition of multipoles of high order in I at the hot spot. For sufficient (dipole) flux to 
be trapped at Rx = 7.5 and 10R* to support a balanced X-wind/funnel flow, our (arbitrary) 
choice of the stellar field thus implies a total field in the hot spot, composed primarily of 
higher order multipoles, that is much larger than the dipole field alone and squeezes the 
funnel into an extremely small spot on the stellar surface. The consequent small covering 
fractions F h and large required fields B h in the hot spots then compare badly with the ob- 
servations. More suitable choices of the multipole distribution, combined with lower values 
of the mass- accretion rate M D in the disk, can do better if relatively large values of Rx / R* 
are indicated (see §4.3). 

Note however that the product FhBh in our numerical calculations is a reasonable num- 
ber: 33, 63, 85 and 105 G respectively, in each of the cases Rx/R* = 2, 5, 7.5 and 10. 
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This suggests that in the context of our calculations, a simple remedy is to choose the actual 
stellar surface to be higher than the fiducial values drawn in Figures 4 and 5. As discussed at 
the end of §3.2, this is equivalent to (1) changing the adopted surface field configuration, and 
(2) concurrently changing the Rx/R* that a particular solution (specified by $t/3>mx and 
Rk/Rx) corresponds to, while (3) leaving the solution parameters at the X-point (§ t /&mx 
and Rk/Rx) unchanged for the specified position of Rx on the equatorial plane with respect 
to the origin. In particular, a larger surface would both increase F h and decrease B h (see 
end of §3.2), rectifying the anomalous values derived for the two largest R x /R* cases above. 

To illustrate this, consider increasing the stellar radius by a factor of 10/7 in the Rx 
= 10i?* case. While the solution parameters at the X-point remain unchanged from the 
R x = 10R* values (<E> t /<E> mx =1.55 and R k /R x — 0.70), the flow now lands at a co-latitude w 
31.8°, and the X-point now effectively corresponds to Rx = (10x7/10)R* = 7R* (Fig. 4). 
The new B r on the stellar surface and B e in the equatorial plane are shown in Fig. 2. We 
see that the magnetic field on the new surface (in arbitrarily scaled, non-dimensional units) 
is lower by a factor of ~20 compared to its original value; the contribution of the highest 
order multipoles has decreased significantly as well, with the quadrupole (/=3) and octupole 
(1=5) components now dominating. The corresponding hot spot is noticeably larger (Fig. 
4), with Fh = 0.91%. Adopting the same B norm (i.e., unchanged dimensional stellar mass, 
radius and accretion rate), (3 and / as above, equation (7) for Rx = 7/?* with the new Fh 
then implies Bh = 8.8 kG, much smaller than the 91 kG derived from our original surface 
field configuration for a similar R x j R* = 7.5. The new F h and B h are perfectly compatible 
with observed CTTs; indeed, they correspond quite closely to specific CTTs configurations, 
as discussed below. 



4.3. The Specific Cases of V2129 Oph and BP Tau 

Recently, Donati et al. (2007) have used spectropolarimetry techniques followed over 
the rotation cycle of the modestly accreting T Tauri star V2129 Oph to map the 3-D mag- 
netic field on its surface. They find that the funnel flow lands on the star in a complex, 
nonaxisymmetric, pattern of hot spots that is distinctly non-dipolar in its geometry. To 
minimize the derived stellar flux, they assume that the unseen hemisphere of the star is 
antisymmetric in its field distribution to the seen hemisphere. They then reconstruct the 
multipolar field distribution by including spherical-harmonic /-values up to I = 9, and find 
that the dominant contribution is from the the quadrupole and octupole components. The 
observed dipole component corresponds to a field B = 350 G at the magnetic pole of the 
star, whose radius they estimate to be i?* = 2.4 R Q . Hence the star has a putative dipole 
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moment fju = BRl/2 = 8.13 x 10 35 G cm 3 . 

If we adopt their data set, M* = 1.35M , M D = 1O~ 8 M yr _1 , and simply use the dipole 
formula for the magnetopause standoff distance (fT2"|) with the OS95 estimate a t = 0.923, we 
get Rx = 2.67-R*, which is appreciably smaller than the corotation radius of R co = 6.67-R* 
corresponding to a 6.53 day rotation period. Part of the difficulty may lie in too large an 
assumed accretion rate for this fairly evolved pre-main-sequence star. From a calibration of 
Paschen and Brackett line luminosities, Natta et al. (2006) obtain = 5.5 x 1O~ 9 M yr _1 
for this object, consistent with the 4 x 1O~ 9 M yr _1 Donati et al. also infer (but use as 
a lower limit) from the observed Ca II flux of V2129 Oph and the M^-Ca II flux relation 
derived for CTTS by Mohanty et al. (2005). However, since we advocate a disk accretion 
rate that is (1 — /) _1 = 3/2 times larger than the stellar accretion rate, the Natta et al. 
number translates to an that does not differ much from the assumed Mr, = 1O _8 M 
yr _1 , and we shall continue to use this value as a standard of comparison. 

With their input data, Donati et al. (2007) claim that the pure (aligned) dipole model of 
Cameron & Campbell (1993) gives values of R m equal to 5.1-R* (magnetic buoyancy model) 
to 8.1-R* (turbulent model), in better formal agreement with the empirical value of R co than 
the pure (aligned) dipole models of Konigl (1991) or OS95. This claim is accurate; however, 
we note that the superior agreement results by choosing a value for the slip parameter 7 
equal to unity, whereas footnote 1 explains that 7 is probably small compared to unity if we 
accept the estimates of Cameron & Campbell at face value. 

The difference between pure dipole and more complex multipole models provides a better 
resolution. Suppose that we can apply the axisymmetric equations of §2.2 to V2129 Oph. 
For M* = 1.35M , M D = 1 x 1O" 8 M yr" 1 , and R* = 2.4R©, equation flU yields B noim = 
27.2 G. With (3 = 1.21, / = 1/3, and R x /R* = 6.67, equation © yields F h B h = 79 G, in 
reasonable agreement with the observed value F^B^ pa 100 G (where pa 5% and B^ w 2 
kG is the mean strength in the hot spots of the total field, not just the dipole component). 

Donati et al. (2008) have carried out similar observations of the CTTS BP Tau, with 
parameters M* = O.8M , R* = 1.95i? , M D w 3 x 1O~ 8 M yr _1 , and rotation period P = 
7.6 days corresponding to Rx/R* = 7.5. They find a strong, mainly axisymmteric poloidal 
field, with roughly equal mean field strengths of 1.2 kG and 1.6 kG respectively in a dipole 
and an octupole component (both slightly tilted with respect to the rotation axis). From 
the stellar mass, radius and accretion rate, equation (8) gives B norm = 53.7 G. With adopted 
(3 = 1.21 and / = 1/3 as before, equation (7) then predicts F h B h = 170 G, in excellent 
agreement with the data: Donati et al. measure a mean hot spot field strength B^ pa 9 kG 
(including both dipole and octupole components) and covering fraction Fh ~ 2%, leading to 
an observed F h B h pa 180 G. 
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It is also noteworthy that the V2129 Oph and BP Tau observations resemble our scaled 
calculation presented at the end of the last section (§4.2). In both the model and the data, 
the surface fields are dominated by multipole components up to an octupole, the Rx are 
similar (6.7-7.5/2* in the data, intermediate 7R* in the model), and the flow lands at high 
latitudes; under the circumstances, it is heartening that (with adopted i?*, M* and Mr, 
comparable to the data) the predicted Fh (~1%) and Bh (8.8 kG) are also quite similar to 
the observed values (2-5% and 2-9 kG respectively). While we have made no attempt to 
exactly reproduce the data (and indeed cannot, given our assumptions of axisymmtery and 
aligned field), it is thus evident that both our analytic mutipole equations and corresponding 
numerical models can reproduce the broad qualitative and quantitative features of the data. 

Thus, the time has probably come to abandon formulae for disk truncation radii based on 
simple dipole models for the unperturbed stellar field in YSOs. We take encouragement that 
progress in both the observational and theoretical developments allow researchers to argue 
meaningfully about missing factors of 2 or so in as complex a subject as the magnetospheric- 
disk interactions of CTTSs. 

On the other hand, the 3-D simulations of Long et al. (2007, 2008) always disrupt the 
disks at a point where "ram pressure balance" holds approximately (see also the calculations 
of R03, R04, Long et al. 2005). However, as discussed in §1, none of these simulations 
produce an X-configuration of trapped flux. Trapped flux creates fields near the inner edge of 
the disk that depart appreciably from an unperturbed stellar dipole. Thus, until simulations 
of the type performed by R08 are carried out with complex (non-dipolar) fields on the stellar 
surface, we regard factors of 2 uncertainties in the location of the stellar magnetopause as 
an open question. 

As a final comment, we should note that although the global field geometries recon- 
structed as Figs. 15-16 of Donati et al. (2007) and Fig. 15 of Donati et al. (2008) are 
visually compelling, the extrapolation of exterior magnetic fields from measured ones on a 
stellar surface is, in general, a mathematically ill-posed procedure and does not produce 
unique configurations. As an example, the boundary condition (16a) of this paper does not 
suffice, by itself, to fix the solution of the elliptic PDE ffTBI) for r > R*. Also needed are 
the other conditions (16b)-(16g). Other models for the disk- magnetosp here interaction will 
involve different outer boundary conditions and/or singular surfaces, lines, and points. The 
combination of empirical measurements of the 3-D fields on the stellar surface, accurate the- 
oretical modeling and better determinations of accretion rates is needed to produce further 
progress in the field. 



5. Conclusions 



Since its inception more than a decade ago, generalized X-wind theory has made many 
successful predictions and survived several key confrontations with observations (e.g., Johns 
& Basri 1995; JG02, Shang et al. 2004; McKeegan 2006; Pyo et al. 2006; Zolensky et 
al. 2006; Carr 2008; Edwards 2008). Claims of large rotation rates seen in YSO jets have 
been used to argue against X-winds (e.g., Bacciotti et al. 2002, Coffey et al. 2004). These 
claims turn out to have difficulties, e.g., indications in some cases that the inferred rotation 
in the jet is counter to that of the disk, or the absence of any obvious rotation in jets best 
oriented (i.e., in the plane of the sky) to show such rotation (e.g., Cabrit et al. 2006, Pety 
et al. 2006, Coffey et al. 2007, Lee et al. 2006, 2007). Present upper limits on the rotation 
in the latter cases show that if any disk wind is present in the observations, they must 
have small launch radii that make them look very similar to X-winds (see Cai et al. 2008 
for a fuller discussion). In a similar manner, the finding that the magnetic fields on the 
surfaces of CTTSs are locally strong as seen in hot spots, but globally weak when averaged 
for the photospheric polarization signal (Valenti & Johns- Krull 2004), turns out not to be 
an argument against disk- locking via funnel flows (Matt & Pudritz 2005), but evidence for 
the phenomenon of flux trapping predicted by X-wind theory (OS95) and quantified in this 
paper (see also JG02). 

In hindsight, we should perhaps not be surprised that fully convective CTTSs have 
surface distributions of magnetic field that are non-dipolar. Such configurations complicate 
the theoretical modeling of funnel flows, but they break the degeneracy in the very differ- 
ent explanations offered by Ghosh & Lamb (1978, 1979a, b; see also Konigl 1991) and by 
Cameron & Campbell (1993), Shu et al. (1994a), OS95, and this paper for the standoff 
distance Rx of the stellar magnetopause where the disk is truncated at its inner edge and 
to whose Keplerian rotation the star is locked in steady state. In particular, this paper 
shows that the observed field strengths in funnel-flow hot spots and their fractional covering 
fractions F h have a rational explanation in the context of generalized X-wind theory, while 
they are inexplicable in other semi-analytic formulations of the steady-state problem. The 
strong concentrations of magnetic flux available to realistic configurations of surface fields 
on young stars may mean that it is harder to crush their magnetopsheres than has been 
estimated by naive dipole estimates. 

While the results from numerical simulations (e.g., Goodson et al. 1999, Kiiker et al. 
2003, Long et al. 2005) contain many of the elements seen in the semi- analytic theory, 
complete agreement has not yet been achieved, perhaps because the numerical simulations 
invariably contain, so far, too much resistive diffusion of the field, relative to angular momen- 
tum transport by turbulent viscosity, to give a good semblance of the phenomenon of trapped 
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flux (for progress in this regard, see Romonova et al. 2008, p. 281). Nevertheless, the trends 
are promising, and at some stage, the awesome computing power of modern machines will be 
able to reliably extend the solutions to the non-axisymmetric and time-dependent regimes 
occupied by actual systems that are inaccessible to semi-analytic techniques. 

We thank Suzan Edwards for informative conversations, Hsien Shang for providing Fig- 
ure 1 and calculations of the wind interface, and the anonymous referee for detailed comments 
that greatly improved the paper. SM acknowledges support from a Spitzer Fellowship; FS, 
from a NSC grant awarded to TIARA in Taiwan, where this work was completed but not 
written up until both authors were back in the USA. 

A. Solution of Governing PDE with Associated Boundary Conditions 

We wish to solve the elliptic PDE (TT5T) subject to the boundary conditions (16a)-(16g). 
Condition (16a) requires that the field approach a sum of aligned multipoles near the origin. 
The goal of our analysis is to mimic the observed characteristics of T Tauri stellar surface 
fields (§3.2). The required sum is thus found by specifying the unperturbed multipole field 
at a finite radius from the origin. A non-restrictive multipole sum has the advantage 
that it defines an analytically differentiable field over all space excluding the origin (see 
§3.2). This allows us to perform the numerical computation without the added complexity 
of introducing the curved stellar surface as a boundary. The actual location of the stellar 
surface is important, however, for astrophysical interpretations (see §3.2). 

Condition (16b) demands that the field take the shape of a complete fan near the X- 
point, where a specified amount of flux $ t is trapped. OS95 show that a fan is the simplest 
curl-free solution for trapping flux at a single point. 

Condition (16c) says that the fields under consideration are confined by the X-wind, 
which cannot be approximated accurately as a vacuum field. Thus, the uppermost wind 
streamline corresponds to the flux tube just beyond $=$ u ,. All the field lines participating 
in the wind (as well as in the funnel flow) have their disk footpoints at the X-point, where 
the field assumes the shape of a fan. Moreover, for a cold flow, the wind streamlines must 
exit the X-point at a shallower angle to the disk than the critical equipotential, in order for 
outflowing material to be successfully loaded onto them. The critical equipotential on the 
wind side of the flow is inclined at $ = 60° to the disk at the X-point. Thus, we take the 
uppermost wind streamline z w (w) to attach to the X-point at an angle of exactly 60° to the 
disk. These considerations enable us to identify the constant flux along z w (w) as § w = $f 
= $i^/7r, where $ w = ir/3. 
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H. Shang (2005, private communication) kindly provided us the locus of the wind in- 
terface z w {w) that she obtained from a detailed numerical calculation of the force balance 
across the dead- zone/X- wind interface, as elucidated by Shu et al. (1995). This locus has 
a slope = tan[7r/3] (= tan^ w ) near the X-point, and approaches a vertical orientation more 
slowly than OS95's more arbitrary parametrization. We fit a quadratic polynomial to the 
numerically obtained locus by z w (zu) = a {w / Rx) 2 + b (zu/Rx) + c where a = 0.197 and b = 
1.333, with c = — (a + b) guaranteeing that the streamline intercepts the disk, z—0, at the 
X-point w = Rx- 

Condition (16d) arises from the following considerations. Since OS95 adopted an ana- 
lytic expression for z w {w) that can be extended to infinity, they were free to choose an upper 
boundary to their computational region at z nmx = oo, and assume there the simple boundary 
condition dQ/dz = 0. In our case, the locus z w (zu) from Shang is numerically calculated and 
thus forced to be spatially finite. Within our computational domain - the region interior 
to the wind interface - we therefore replace the OS95 condition of a uniform distribution of 
longitudinal field at z max = oo by one at a large but finite z max = 15. The flux function that 
gives a uniform value for B z is $ = Q w (w /w niax ) 2 , where we have imposed the boundary 
values that $ = on the axis w = and $ = $ w = $ t /3 at w = ro max when we reach the 
wind interface at the height z = z max . 

If all the disk material came to an abrupt end at the X-point, then field lines fanning out 
from the X-point down to $ = n (i.e., down to the midplane interior to the X-point) would 
reconnect with their oppositely directed counterparts across the midplane. To prevent this, 
we (like OS95) relax the constraints to allow a current sheet (later called "the reconnection 
ring") to extend slightly inward of Rx, to the radius Rk- The latter, the inner terminal 
point of a current sheet (which extends in its exterior into the disk for w > Rx), is called 
the "kink-point" radius and is the vertex of a magnetic "Y" if we consider the corresponding 
field geometry below the mid-plane (see below). The last horizontal fan flux-tube from 
the X-region (both above and below the midplane), with value $ = $ t and carrying the 
terminal part of the funnel flow, is then able to skim the midplane between Rx and Rk 
before bifurcating into two lines that form the innermost boundaries of the out-of-plane 
funnel flow. This bifurcation encapsulates the statement contained in condition (16e). 

Inward of Rk, the inspiraling material is forced out of the disk plane. The magnetic lines 
of force encountered in the midplane at these radii are closed multipole fields with footpoints 
in the star but not the disk. For decreasing radii, they become increasingly unperturbed as 
the extension of the vacuum fields above the star. Reflection symmetry across the midplane 
then gives rise to condition (16f). Note that conditions (16e) and (16f) imply that the field 
at the midplane is completely vertical for w < Rk , and completely horizontal for w > Rk- 
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The field at R k must then have a kink, as illustrated in Figure 1. We see that $ = $ t thus 
represents the last field line emanating from the star that is pinched into the X-region, with 
a long slender beak extending from R k to Rx- OS95 demonstrate that near R k , the ideal 
geometry of the magnetic field lines in the meridional plane is an even three-way split of the 
full angle 2ir, among the three magnetic sectors - field lines leading towards the X-point, 
emanating from the X-point, and unconnected to the disk - that have equal field strengths 
across their interfaces (Fig. 1). Thus, the last flux line connected to the X-point, $ t , subtends 
an angle of 120° as it leaves the midplane at R k (creating the kink). In our (and OS95's) 
formulation, the kink radius R k represents an eigenvalue of the problem, whose value must 
be determined as part of the solution for B p . In other words, R k must be specified such that 
the derived B p conforms to the physically valid field geometry near R k described above. 

Finally, boundary condition (16g) sets the zero of the flux- function along the z-axis for 
axial symmetry. Thus, the field line exiting the pole of the star ($ = 0) attaches to its open 
counterpart in the wind region only at z = oo. 

Conditions (16c)-(16g) specify $ or its normal derivative on a closed curve. For a 2-D 
elliptic PDE in $ such as equation (1151) . these alone usually suffice to completely determine 
$ throughout the bounded region. Yet we (and OS95) are able to impose, apart from 
the midplane Y (which has an associated eigenvalue), two further conditions, (16a) and 
(16b), at the origin and X-point respectively, without over-determining the problem. This 
is because the required solution for $ is not differentiable at the origin (singular dipole 
for OS95, singular multipoles for us) or at the X-point (magnetic fan). The multipole 
and fan behaviours must be imposed explicitly in order to obtain the desired solution that 
implicitly has current sheets and divergent current loops at the boundaries and the origin 
of the problem. To deal with these singularities at the origin, X-point, and midplane Y, we 
follow OS95 in solving the PDE as follows. 



where $ m and $f are the already defined multipole and fan flux functions; $ r is the remainder 
field that must be computed to completely determine the final solution $. We also define 
the differential (Stokes) operator, 



so that the PDE of equation (fT5l) is equivalent to £>($) = 0. Since the multipole field $ 
itself obeys T>($ m ) = 0, the remainder field must, by equation flAll) . satisfy 



We define 



(Al) 




(A2) 



T?($ r ) = -P($ f ). 



(A3) 
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Combining the earlier boundary conditions on $ (eqs. 16a-g) with equation flAll) . we derive 
the corresponding boundary conditions on $ r to be 



$ r = — $ m on z w (w) 

$ r = $ f (ro/u7 max ) 2 - $ f - $ m onz = z max 
$ r = -$ m on z = for R k < vo < R x 
d^/dz = -d${/dz on z = + for < w < R k 
<I> r = -$ f on = 



(A4a) 
(A4b) 
(A4c) 
(A4d) 
(A4e) 



We are interested in the geometry of the total field. Because the governing PDE is linear 
and equidimensional, the geometry is unaffected by spatial or field-strength scaling. All 
spatial parameters are thus normalized by i?*, and field-strengths by $ mx , the unperturbed 
multipole flux at the X-point in the absence of interactions with the disk. That is, the 
problem is scaled so that = 1 and $ mx = 1. With the above boundary conditions, 
solving equation ( 1A3I) for the scaled $ r then requires only a specification of two quantities: 
(i) $ 4 /$ mx , the ratio of flux trapped at the X-point to the unperturbed multipole flux there 
in the absence of a disk; and (ii) Rk/Rx, the ratio of the kink point radius to the X-point 
radius. The value of $t/$ mx is determined as the minimum required to have no "frustrated" 
field lines for steady funnel flow (i.e., no field line leaves the X-point on the "downhill" side of 
the critical equipotential, yet subsequently crosses it again so that a cold funnel flow would 
encounter an insurmountable potential barrier; see OS95), while the eigenvalue R k /Rx is 
found by trial and error, given $t/$ mx , by requiring the solution to satisfy the ideal Y- field 
geometry near R k . The derived remainder field $ r is then superimposed on the (scaled) $ m 
and $f to obtain the total scaled solution $ from equation flAll) . 

To solve equation flA3p and its associated boundary conditions on a grid, we use the 
transformation {vj,z) — > {w, C), where ( = 2\[yJ(z/\) - 1-1]. Here A = (2a + 6) 2 /4a, where 
a and b are the coefficients of our quadratic form for the wind interface defined earlier. This 
transformation maps the original curved wind interface z w {w) = a w 2 + bw + (a — 6) to the 
inclined straight line = 2\fa\[w — 1], yielding a trapezoidal grid. The PDE is then 

solved by simultaneous overrelaxation through finite-differencing. 
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Fig. 1. — Original X-wind accretion/outflow solution, using a stellar dipole field, that com- 
bines the funnel flow of OS95 with the X-wind solution of Shu et al. (1995) for the case Rx = 
5.R* (modified version of figure kindly supplied by H. Shang). Solid curves are magnetic field 
lines: accretion funnel flow lines in red, open field lines from the X-point in blue, all others in 
black. Key features generic for good solutions: flux-trapping at the X-point Rx, extension of 
an infalling disk slightly inward of Rx until the kink-point R& which makes an angle of 120°, 
and the uppermost funnel flow field line always remaining below the equipotential curve. 
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But notice the high net field polarization on the stellar surface, and the relatively large hot 
spot covering fraction (~6% if R x /R* = 5) implied by the dipole field. See §§2 and 3.1 
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Fig. 2. — The unperturbed stellar magnetic field (B r and Bg) corresponding to the $ m of 
equation (25), in non-dimensional, arbitrarily scaled units (but same relative scaling in all 
panels). Top left: B r on stellar surface, as a function of fi. The various multipole components 
I = 1-11 are shown as colored lines (see key in plot); the corresponding total B r is in black. 
Top right: Total Bg on stellar surface as a function of [i (individual multipole components 
not shown). Both B r and Bg oscillate rapidly, signifying a loopy field structure on the 
surface. Middle left: Bg as function of radius (in units of stellar radius) in the equatorial 
plane (with R* = 1; B r — in this plane). Individual multipole components in color (same 
key as in top left panel), total in black. Middle right: Zoom-in of total Bg in the equatorial 
plane at large radii. Note that Bg changes direction at l.Oli?* and 3.25-R* (passing from 
negative to positive values), finally decreasing smoothly to zero at large radii as the dipole 
field dominates. Bottom left: B r as a function of /x on rescaled stellar surface, R* = 10/7 (as 
opposed to R* = 1 in the previous panels). Bottom right: Bg as a function of radius (in units 
of stellar radius) in the equatorial plane, for the same rescaled stellar surface. Multipole 
components in both panels in color (same key as in top left panel), total field in black. From 
comparison to top left and middle left panels, note that (1) B r and B g are ~20 times weaker 
for the larger rescaled surface compared to the smaller original one, and (2) the quadrupole 
and octupole (/ = 3 and 5) components largely dominate for the rescaled surface. See §§3.2 
and 4 (particularly §4.2 for the rescaled surface). 
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Fig. 3. — Magnetic lines of force in the meridional plane, corresponding to the $ m of equation 
(25). Note that this $ m is non-dimensional and arbitrarily scaled; absolute scaling and 
dimensional field strengths are discussed in §4.2. Top panel: Zoom-in of field close to the 
star. Contours trace lines of constant flux, from <£> m = -0.035 to +0.035, at intervals of 0.0005. 
The field is highly loopy near the stellar surface. Bottom panel: Same as above, except over 
a larger spatial range. Contours are from -0.002 to +0.002, at intervals of 0.00002 (1/15 of 
the range and 1/25 of the spacing in the top panel). Far from the star, the field increasingly 
resembles a simple dipole. Colors denote contour levels: blue — > red =>- -ve — > +ve. The field 
direction thus corresponds to progression in colours with r or 6 (for Bq and B r respectively). 
Red asterisks mark the two null points in the field in the equatorial plane, where B e changes 
direction. See §§3.2 and 4. 
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R x = 10 R„ Rk/Rx = 0.70, * t /*mx = 1.55 




2 4 6 8 10 12 14 

w/R. 



Fig. 4.— Final field solution for R x = 10R*, with R k /R x = 0.70 and $ t /$ mx = 1.55. 
Contours are lines of constant flux. Funnel flow field lines in green, all others in black. 
Dashed line shows the wind interface. Equipotential curve shown in red: the uppermost 
funnel flow line always lies below this curve, ensuring steady accretion. Angle of 120° to 
the disk at the kink point Rk shown in blue: this is the angle at which the last funnel flow 
line exits the disk at R k - The yellow curve represents the stellar surface i?* = 1, on which 
the original unperturbed field is specified. The field geometry near Rx closely resembles the 
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dipole solution in Fig. 1, but the funnel is squeezed into a much narrower flow near the 
stellar surface, producing a very small hot spot (covering fraction 0.065%) compared to the 
dipole case. The brown curve represents a rescaled stellar surface R* = 10/7; the Rx = 10R* 
solution shown here for the original surface (in yellow) becomes the Rx = 7.R* solution for 
this new surface. The hot spot covering fraction for the new surface is much larger (~1%) 
than for the original one. See §§4.1 and 4.2. 



= 7.5 R„ Rk/Rx = 0.61, * t /*mx = 1-71 




Fig. 5.— Final field solution for R x = 7.5i?*, with R k /Rx = 0.61 and $t/$ mx = 1.71. 
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Colors and symbols same as in Fig. 4. See §§4.1 and 4.2. 

Rx = 5 R„ R k /R x = 0.247, * t /* m x = -25.0 




Fig. 6.— Final field solution for R x = 5R*, with R k /R x = 0.247 and $ t /$ mx = -25.0. 
Colors and symbols same as in Fig. 4. Note that the hot spot has now jumped to a different 
location, with the opposite polarity in B r , compared to Fig. 5. See §§4.1 and 4.2. 
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Fig. 7.— Final field solution for R x = with R k /R x = 0.87 and $ 4 /$ mx = 40.0. Colors 
and symbols same as in Fig. 4. See §§4.1 and 4.2. 
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Table I. Parameters of Unperturbed Stellar Multipole Field 



I A l /A 1 



1 


1 


3 


+ 1.64 


5 


-3.46 


7 


+0.82 


9 


+ 1.45 


11 


-0.72 



